Effect of weak rotation on large-scale circulation cessations in turbulent convection 
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We investigate the effect of weak rotation on the large-scale circulation (LSC) of turbulent 
Rayleigh-Benard convection, using the theory for cessations in a low-dimensional stochastic model 
of the flow previously studied. We determine the cessation frequency of the LSC as a function of 
rotation, and calculate the statistics of the amplitude and azimuthal velocity fluctuations of the 
LSC as a function of the rotation rate for different Rayleigh numbers. Furthermore, we show that 
the tails of the reorientation probability distribution function remain unchanged for rotating sys- 
tems, while the distribution of the LSC amplitude and correspondingly the cessation frequency are 
strongly affected by rotation. Our results are in close agreement with experimental observations. 
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The complex phenomenon of thermal turbulence, also 
known as turbulent Rayleigh-Benard convection (RBC), 
develops in a heated fluid layer under gravity by a succes- 
sion of instabilities in the thermal transport due to an in- 
terplay between different driving forces such as buoyancy, 
viscous drag, and diffusion. This balance between differ- 
ent transport mechanisms is quantified by the Rayleigh 
number Ra — a^^g/STL^ / vk, determining the flow state. 
Here, a^, is the isobaric thermal expansion coefficient, g is 
the gravity field, AT is the temperature gap between bot- 
tom and top plates, L is the height of the fiuid container, 
K is the thermal diffusivity and v is the kinematic viscos- 
ity. When Ra is sufficiently large, the flow becomes tur- 
bulent and a large-scale circulation (LSC) is formed. The 
latter is maintained by the emission of plumes from the 
top and bottom surfaces which, due to buoyancy, move 
upwards (hot plumes) or downwards (cold plumes) [IH3] ■ 
This LSC is known to appear in various rotating natural 
systems, such as atmospheric \^ and oceanic flows |5j, 
and the dynamo driving planetary magnetic fields [6]. 

Such a circulating state does not persist indefinitely, 
however, and cessations followed by restarted flows at a 
new azimuthal angle occur sporadically [U [7HTU]. The 
complex dynamics of the LSC is influenced by the heat 
and momentum transport mechanisms, as well as by the 
geometry and aspect ratio F (diameter over height) of 
the fluid container. In experiments using a cylindrical 
geometry with F = 1, the LSC occurs in a nearly vertical 
plane, whereas for aspect ratio F — 0.5, two convection 
rolls may coexist [TT] . Recent studies also report that the 
LSC flow-reversals and its tilted orientation are strongly 
influenced by the corner-flows that form in a rectangu- 
lar geometry |12j . The dynamics of a nearly- vertical, 
single LSC can be modeled by a set of nonlinear stochas- 
tic differential equations that describe the amplitude of 
azimuthal temperature variations, 5, and the azimuthal 
orientation angle, [131 HI- This model is found to be 
in excellent agreement with the typical fluctuations of the 
LSC J^il4j. Yet, the cessations require an extension of 
the model, since boundary terms describing the thermal 
and viscous diffusion become dominant terms when the 



amplitudes are small, as they inevitably must be during a 
cessation. With this extension, the model is able to pro- 
vide an excellent description of the cessation rates, as well 
as the LSC amplitude and azimuthal velocity probability 
distribution functions (PDFs) including the tails [15) . 

The purpose of this Letter is to study the effect of 
weak rotation on the statistics of the main degrees of 
freedom, i.e. amplitude and azimuthal velocity, and the 
LSC cessation events in RBC. Here the rotation rate Vl is 
conveniently represented by the dimensionless convective 
Rossby number, Ro = {2n)~'^ a^g/^T / L, which mea- 
sures the buoyancy relative to the Coriolis force. The 
total heat transport, measured by the Nusselt number M 
relative to heat diffusion, has a non-trivial dependence 
on Ro, exhibiting three different regimes [101 [HI [T7] : (i) 
an LSC-dominant regime at Ro"^ < Ro;:^ ~ 0.4 where 
buoyancy dominates over Coriolis effects, (ii) an Ekman- 
vortex pumping regime at Ro^^ < Ro~^ < Ro"^^. where 
the Coriolis force dominates over the buoyancy force for 
a non- vanishing Ro~^, and (iii) at Ro~^ > Ro~^^ (typ- 
ically Ro~^j. ~ 3 — 10 depending on Ra [10]) where the 
velocity gradients, and thus heat transport by convection, 
are suppressed via the Taylor-Proudman effect [10] . Here 
Ro""'^ scales as l/F due to a finite-size effect [El [19], while 
Ro~^j. is determined by an interplay between the Ekman 
enhancement and the Taylor-Proudman depression. 

In this work we focus only on weak rotations, where the 
LSC is still present, namely in the regime < Ro^^ < 0.6. 
Recent experimental studies report a non-monotonic de- 
pendence of the LSC mean amplitude on Ro~^ |10j . 
The mean amplitude initially increases with increasing 
Ro^^ until Ro^^ becomes close to (but below) a criti- 
cal value, Ro^^. The origin of this increase is not well 
understood, but may be associated with centripetal ef- 
fects [20], see below. For higher Ro~^, the detaching 
plumes from the thermal boundary layers (BLs) interact 
with the Ekman- vortex structure that forms at RoJ^, so 
that these plumes can no longer feed the LSC jH]. Thus, 
the mean amplitude decreases at Ro^^ > Ro~^. For 
Ro^^ > 0.6 — 1 (depending on Ra, Pr and F), the mean 
amplitude becomes comparable to the root-mean-square 
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temperature fluctuations about the mean, and thus be- 
comes iU-defined [10]. This non-monotonic behavior of 
the LSC mean ampHtude on Ro"^ is accompanied by a 
non-monotonic behavior of the LSC cessation frequency, 
and this is what we calculate below. This dependence 
contrasts with our result for the tails of the reorientation 
statistics, namely the probability for a large angle change 
in the LSC plane, which we show is almost independent 
of the rotation strength in this weak rotation regime. 

Statistics of the LSC dynamics with rotation:- The evo- 
lution of the dimensionless LSC amplitude, ^, and az- 
imuthal angle, 9o, in the absence of rotations is governed 
by the stochastic differential equations [15] : 
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Here, ^ is a constant related to the heat transport or 
the inverse of the thermal BL width, ^ — S/Sq is the di- 
mensionless LSC amplitude, where S is the physical LSC 
amplitude, Sq « ATcrRe'^^^/Ra is the mean LSC ampli- 
tude, and a = v/k is the Prandtl number. Furthermore, 
the Reynolds number is defined as Re — (ra/rg)^ :s> 1 
where ts and Tg are the turnover times in the LSC and 
azimuthal planes, respectively, while time in Eqs. ([T]) is 
measured in units of the corresponding turnover times ts 
and Tg [15]. Also, delta-correlated Gaussian stochastic 
forcing terms fs (t) and fg (t) with amplitudes Dg and Dg 
are included in Eqs. (jlj to simulate the effect of turbulent 
fluctuations. Finally, the coefficients a,/3, = 0(1) 

are included to account for the geometric coefficients 
from the spatial volume averaging procedure |131 I14j . 
Note, that in Eqs. (fTl) all tildes (which appeared in [TS] 
due to time rescaling) were removed for clarity. 

In the presence of rotation, the equation for the LSC 
amplitude [the first of Eqs. ([T])] remains the same, since 
the same drag and bouyancy forces drive the motion of 
the LSC in the vertical plane. Yet, the coefficients in 
front of these terms may depend on Ro, to account for 
the experimental fact that the mean dimensionless LSC 
amplitude is a function of Ro~^. To find this Ro- 
dependence, we notice that the coefficients Q!(Ro), /3(Ro) 
and A(Ro) are related to each other by the constraints 
that the PDF P(^) is centered around (, — £,o and has a 
width Ds 15]. Employing these constraints, we obtain 

a(Ro) = 1 - 3A/eo , /3(Ro) = ^0-^/' (1 - 2A/eo) . (2) 

Now we formulate a simple theory that can explain the 
dependency of these coefficients as well as that of on 
51. We assume that a = uq, is constant (because the 
geometrical coefficient in the buoyancy should not be af- 
fected by rotation) and expand /3 — /3o(l + &ifl + &2^^^) 
(assuming that the thickness of the viscous BL is de- 
pendent on 51), where ao = 1 — 3^0i = 1 ^ SAq 
and Aq — A{fl — 0). In this way, we can find ex- 
pressions for A{Q) and Co(^^)i by equating these expres- 
sions for a and (3 with Eq. ([2|, and keeping terms up to 
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FIG. 1. (Color online). Shown are the values of (triangles), 
and A (squares). The dashed line is a parabolic fit. The up- 
per inset shows /3 (triangles) compared with (squares), 
while the dashed line is a parabolic fit. The lower inset shows 
a and the dashed line is a guide for the eye. All quantities 
are normalized by their non-rotating values. 
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and A{il) = ^o'?o(51)- Note, that plugging this into 
Eq. we obtain a(Ro) — ao and /3(Ro) ~ Po^.o'''^^- 
Because /3{Ro) represents the inverse Rossby-dependent 
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width of the viscous BL, the latter increases as for 
< Ro"^ < Ror\ in agreement with recent experimen- 
tal observations [Ml jH] . The origin of the increase in the 
viscous BL width in the weak-rotation regime can be as- 
sociated with centripetal effects [20] that tend to increase 
the BL width according to the Prandtl-Blasius theory. 

In Fig. [T] we plot the experimental mean amplitude 
and A (normalized by their value at zero rotation) as a 
function of Ro^^. According to our theory, provided that 
^0 is well- approximated by a parabolic expansion in 51, 
in the weak rotation regime, A/Aq should coincide with 
Co/Co(0), and this is indeed the case. The insets in this 
figure confirm the dependencies of a and /3 on 51 for two 
different Ra numbers. Note that in order to extract A 
from the experimental data, we have used the relation 
A = BDs/2 [15] where B and Ds are specified below. 

Since the equation for the LSC amplitude [the first of 
Eqs. ([!])] is independent of the azimuthal reorientation 
angle 6o, we can analyse it separately by writing the cor- 
responding Fokker-Planck equation and finding its sta- 
tionary solution 15J. Defining the potential 
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the frequency-dependent PDF, P{£,, Ro) takes the form 
P(e,Ro) = Ce-2^(«)/^^ (4) 

Here C is a normalization constant, Ds = Ds(Ro) is 
the PDF width in the Gaussian regime: w ^o) ^ 

g-(C-4o) /i2Ds)^ g^jj^j Q _ 5(j|q) is the logarithmic deriva- 
tive of the PDF at small ^, since < 1) e^«. 
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FIG. 2. (Color online). In panels (a)-(f) shown are PDFs of 
the LSC rescaled amplitude ^ for different Ro and Ra num- 
bers. Each row represents a different Ra number. The exper- 
imental PDFs are represented by triangles, while the fitting 
curves (dashed line) are the analytical PDFs Q with param- 
eters B and Dg determined experimentally for each Ro num- 
ber. The width of the PDFs initially increases with increasing 
Ro~^ while the slope of the left tail decreases. This can be 
seen in panels (g) and (h) where we plot the experimental B 
and Ds as functions of Ro~^ for different Ra. 



Furthermore, it has been shown that the cessation fre- 
quency - the frequency of events that the LSC amplitude 
goes below a threshold amplitude ^min ^ 1 - is given by 
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-2Z5-i[y(c*)-v(i)]^ 
. (5) 

Here T{^*) is the mean time it takes the amplitude to 
reach ^* ^ 1. 

In order to compare the results for the PDF [Eq. Q] 
and cessation frequency [Eq. ([5])] with experimental re- 
sults, we extract the values of B{Ro) and Ds{Ro) from 
the experimental PDFs, just as was done for the non- 
rotating case [15j . In Fig. [2] we compare experimental 
and theoretical PDFs for different Ro and Ra numbers, 
in the weakly-rotating regime. The theoretical predic- 
tions hold well for various rotation frequencies, provided 
that we use the corresponding frequency-dependent pa- 
rameters. As Ro~^ ~ is increased, the width of the 
PDF increases compared to the non-rotating case while 
the slope of the left tail decreases. This functional depen- 
dence of the experimental B and Dg on Ro^^ is shown 
in the lower panels of Fig. [2] for two different Ra number. 

Finally, in Fig. [3] we compare the theoretical and ex- 
perimental results for the cessation frequency as function 
of Ro for different Ra values and observe good agreement. 
Here, the threshold for cessation was chosen to be 0.15. 

Rotation rate of the azimuthal plane:- In this section we 
investigate the effect of weak rotation on the PDF for 
the angular velocity of the azimuthal plane [the second 
of Eqs. ([ij]. Rotation brings about a Coriolis force which 
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FIG. 3. (Color online). Cessation frequency (normalized by 
its non-rotating value) as a function of Ro"^ for different val- 
ues of Ra number. The triangles are experimental results, 
while the squares are the theoretical predictions according to 
Eq. |[5|. The cessation threshold chosen here was 0.15. 



is proportional to 51 x V^, where \ Vg \ ~ ^ [T3] . As a result, 
the effect of rotation on the azimuthal dynamics amounts 
to the addition of a term of the form ^17$(6'o) to the 
equation for 6q, where $ = 0(1) is some function of the 
azimuthal angle [52]. Therefore, the second of Eqs. ([T]) 
under rotation becomes 



+m- (6) 



As expected, the addition of the first term on the right- 
hand-side of Eq. ([6| changes the steady state solution for 
at ^ ~ 1 from 9 — a.t zero rotation, to ^ il at 
n > [since $(6lo) = 0(1) and Tg/rg < 1]. This has 
been experimentally observed by Zhong and Ahlers [TU] 
in the weakly-rotating regime, below Ro;:^ The reason 
for this limitation is that the steady-state solution is valid 
only as long as the parameters involved are of order unity, 
which holds for not too large f2, but breaks down for high 
rotation frequencies. In the latter case, when Ekman 
vortices start to form, the model becomes invalid. 

Apart from the steady-state solution that represents 
the mean angular velocity, here we are mainly interested 
in the rare events of large deviations in AO ^ 9. The tails 
of the PDF P(A6') were calculated in [T2] for non-rotating 
systems. Here, we show that these tails are unchanged 
by the rotation. The reason is that the term due to ro- 
tation added to Eq. ^ is proportional to ^. As a result, 
large deviations in A0, which occur when ^ ^ 1, are still 
governed by the term proportional to v^^o in Eq. ([g]), 
which is dominant in the regime of ^ <§; 1 |15j . This term 
does not depend on the rotation frequency. Therefore, 
the tails of P(A6') remain unchanged when rotation is 
introduced, as shown in Fig. [4j 

Indeed, in the non-rotating case the tails have been 
found to scale as a power law with exponent —4, while 
experimental results for non-rotating systems demon- 
strate a power-law tail with exponent ^ —4.3 [TS]. In 
Fig. |4|^b) shown are experimental PDFs, P(A0), for dif- 
ferent Ro numbers. In panel (c) we show PDFs aver- 
aged over a wide range of Ro~^ numbers jlOl (generally 
< Ro^^ < 0.6), for different Ra numbers. The dashed 
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FIG. 4. (Color online). PDF P{Ae) as function of Ae. In (a) 
shown are experimental PDFs for four different Ro numbers 
(see legend in panel b) for Ra = 1.8 x 10^°. As Ro~^ ~ 
Q, increases, the left tail becomes fatter while the right tail 
remains more or less constant. In (b) shown are the same 
PDFs on a log-log plot, demonstrating that the tail is well 
described by a power law independent of Ro; the dashed line 
is a power law with exponent —4.3. In (c) shown are PDFs, 
averaged over a wide range of Ro~^ numbers (generally < 
Ro~^ < 0.6), for different Ra. The dashed line is the same as 
in panel (b). This panel implies that the rare-event statistics 
of P{A6) are universal and independent of Ra and Ro. 



lines are the theoretical prediction: a power law curve 
with exponent —4.3. The excellent agreement between 
theory and experiment indicates that the exponent for ro- 
tating systems remains unchanged compared to the non- 
rotating case. Note, that panels (b) and (c) are shown on 
a log-log scale which only allows us to show the positive 
6 region. We have checked that the left and right tails of 



the PDF scale with the same exponent within 5%. 

Angular- dependent asymmetry of the PDF for azimuthal 
fluctuations:- We conclude by pointing out an interesting 
corollary of our analysis. From Fig. |4][a), it is apparent 
that in contrast to the non-rotating case, the PDF P{A9) 
develops an asymmetric shape as a function of A9 that 
is an increasing function of the rotation rate, up to the 
vicinity of 1/Ro ~ 1/ROc. We speculate that this asym- 
metry in the azimuthal velocity fluctuations is related 
to spiral defects and their preferred motion. Near the 
onset of convection, spiral defects are formed in the pres- 
ence of rotation and their azimuthal motion is against 
the direction of rotation [53] . With increasing Ra num- 
ber, the flow becomes more turbulent and the effect of 
the individual vortices on the azimuthal velocity statis- 
tics diminishes [23] . However, as the strength of the LSC 
decreases when Ro~^ exceeds Ro~^, the flow, hence the 
velocity statistics, is more influenced by the preferred 
motion of the spiral defects in the presence of rotation, 
which would explain the observed broadening of the left- 
tail of the PDF P{A9) corresponding to high fluctuations 
against the direction of rotation. This effect of the left- 
tail broadening with increasing fl can be clearly seen in 
Fig.|4]Ja). However, the mechanism relating the asymme- 
try in P{A6) to spiral defects needs further elaboration, 
and is beyond the scope of this paper. 
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